Transport Measurements of Strongly-Correlated Electrons on Helium in a Classical 

Point-Contact Device 
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We present transport measurements of electrons on the surface of liquid helium in a microchannel 
device in which a constriction may be formed by a split-gate electrode. The surface electron current 
passing through the microchannel first decreases and is then completely suppressed as the split-gate 
voltage is swept negative. The current decreases in a steplike manner, due to changes in the number 
of electrons able to pass simultaneously through the constriction. We investigate the dependence 
of the electron transport on the AC driving voltage and the DC potentials applied to the sample 
electrodes, in order to understand the electrostatic potential profile of the constriction region. Our 
results are in good agreement with a finite element modeling analysis of the device. We demonstrate 
that the threshold of current flow depends not only on the applied potentials but also on the surface 
electron density. The detailed understanding of the characteristics of such a device is an important 
step in the development of mesoscopic experiments with surface electrons on liquid helium. 

PACS numbers: 73.20.-r, 73.23.-b, 85.30.Hi, 45.50.Jf 



I. INTRODUCTION 

An electron close to the free surface of liquid helium 
is subjected to an attractive force due to a weak image 
charge formed in the liquid 1,2 . In the presence of an 
applied electric field E z , the potential in the direction 
perpendicular to the helium surface may be written as 
V(z) = -Ae 2 /z + eE z z where A = (e - l)/4(e + 1), e 
is the dielectric constant of the liquid, e is the absolute 
value of the electronic charge and z is the distance from 
the helium surface. The electron is prevented from en- 
tering the liquid due to a ~ 1 eV potential barrier at the 
surface and so remains localized above the liquid. The 
potential V(z) gives rise to a series of bound states for 
perpendicular motion 3 -; in the ground state, the expecta- 
tion value of z is ~ 11 nm. The lifetime of these bound 
states is predicted to be long, leading to the proposal 
that electrons on the surface of liquid helium may be 
good candidates for quantum bits^. 

For a two-dimensional system of charges on a liquid he- 
lium surface, the electron surface density n s may be var- 
ied over a wide range, up to a theoretical limit given by 
the hydrodynamic instability of the bulk liquid surface 5 - - — 
of ^2.2xl0 9 cm~ 2 . By varying temperature the scat- 
tering processes which determine the electron mobility 
in the plane parallel to the helium surface may be con- 
trolled. Below ~ 0.8 K the density of gas atoms above the 
liquid surface becomes effectively zero and scattering oc- 
curs only with excitations of the liquid surface, ripplons, 
leading to high mobilities^ in excess of 10 8 cm 2 /V-s. As 
the Coulomb interaction between electrons is essentially 
unscreened, and the electron separation is much larger 
than the thermal electron wavelength, SSE have been 
used to study classical effects in strongly-interacting elec- 
tron systems, such as the transition from an electron 
liquid to a 2D Wigner solid as the temperature of the 



electron system is decreased 9 . 

Quasi-one dimensional SSE systems^ have been inves- 
tigated using samples in which liquid helium was confined 
in grooves on a dielectric substrate^. The width of the 
grooves was as low as 1.25 /im. In such experiments, the 
temperature dependence of the electron mobility 1 - 2 - was 
found to be in good agreement with that predicted by 
theoretical calculations where transitions between quan- 
tized energy levels for the lateral motion of the electrons 
were taken into account^. However, for such dielectric 
substrates the mobility was also dependent on substrate 
defects in regions where the helium film was thin. 

Recent experiments have made use of microchannel 
devices fabricated by lithographic techniques in order 
to study the properties of SSE in confined geometries. 
A schematic representation of a set of such microchan- 
nels is shown in Fig. [TJ Placed a distance h above 
the bulk surface of superfluid helium, the channel will 
fill by capillary action; the helium surface may then be 
charged. The radius of curvature R of the liquid is given 
by R — a/ pgh + n^e 2 /2eeq where a and p are the surface 
tension coefficient and density of liquid 4 He respectively, 
g is the acceleration due to gravity and £q is the dielec- 
tric constant of vacuum^. The transport of electrons 
in such devices was first demonstrated for a channel of 
width 30 pm and depth 1 /im 15 . The non-linear trans- 
port of the Wigner solid on helium surfaces in microchan- 
nels of width 8-20 /im was studied in three-terminal de- 
vices comprising of source, drain and gate electrodes 
submerged beneath the helium surfac o 16 ' 17 . The ultra- 
efficient transfer of a small number of electrons along par- 
allel microchannels of width 10 /im, again using a series 
of gate electrodes beneath the helium surface, has also 
been demonstrated^ 8 -. In still more advanced devices, 
small ensembles of electrons, including a single electron, 
trapped in a microfabricated circular pool of radius 10 
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FIG. 1. Schematic picture of a set of microchannels. The 
microchannels are filled by the capillary action of superfluid 
4 He, the surface of which is then charged with electrons. 



fim, have been studied using a charge-sensitive super- 
conducting Single Electron Transistor (SET) positioned 
beneath the helium surface^. A Field Effect Transis- 
tor (FET) for electrons on a thin helium film, where the 
electron density may exceed the hydrodynamic limit, has 
also been demonstrated 2 ^ There the separation between 
the split-gate electrodes was 200 /jm. 

Such experiments demonstrate a progression towards 
the sensitive control and measurement of small numbers 
of SSE on helium, raising the possibility of studying novel 
phenomena associated with classically interacting charge 
systems in confined geometries. Many theoretical stud- 
ies have been conducted on the behavior of such sys- 
tems. In quasi-one dimensional systems, at sufficiently 
low temperatures, particles are predicted to form a series 
of rows, the number of which changes with the parti- 
cle density or confinement strength, leading to structural 
phase transitions and re-entrant melting processe a 21 ' 22 . 
Similar phenomena have also been predicted in circularly 
symmetric parabolic confinements^. The pinning and 
depinning dynamics of charged particles at potential con- 
strictions have also been investigated using Monte Carlo 
calculations 2 ^—. However, experimental difficulties have 
restricted progress towards investigating these phenom- 
ena. SET charge measurements, whilst extremely sen- 
sitive, may be plagued by intrinsic two-level fluctuator 
charge nois o 27 i 28 . Also, the electrostatic potential profile 
in microchannel devices may be distorted by contact po- 
tentials and surface charging effects which can be difficult 
to quantify^. 

As a step towards overcoming these experimental dif- 
ficulties, we have performed transport measurements of 
SSE on superfluid 4 He in a microfabricated device. In 
this sample, two microchannel SSE reservoirs are sep- 
arated by a split-gate electrode which, at appropriate 
bias, forms a constriction. As the split-gate voltage V g t 
is swept negative the current / flowing through the con- 
striction is reduced and then reaches zero at a thresh- 
old voltage V*]?. The appearance of step-like decreases 
in current in this device, where each step corresponds 
to a change in the number of electrons able to pass si- 
multaneously through the constriction, has already been 
reported^. Here we present further transport measure- 
ments, investigating in particular the response of the sys- 
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FIG. 2. Point-contact device for electrons on helium, (a) 
Two arrays of microchannels form the left and right electron 
reservoirs, (b) A split-gate electrode is positioned at the base 
of the small central channel which links the two reservoirs, 
(c) Diagram showing the dimensions of the electrodes in the 
central channel. 



tern to increasing AC driving voltage, which causes sig- 
nificant non-linear transport effects to emerge. These 
results are discussed in relation to an electrostatic model 
of the sample which shows that offsets in electrode po- 
tentials, arising for reasons which are not yet clear, are 
important in determining the potential profile of the de- 
vice. 



II. EXPERIMENTAL 

The sample used in this experiment was prepared us- 
ing multilayer optical and electron-beam lithography on 
a Si wafer, the surface of which was oxidized. Two ar- 
rays of microchannels, defined by a guard electrode, act 
as electron reservoirs between which electrons may be ex- 
changed. These left and right reservoirs consist of 25 mi- 
crochannels of width w = 20 /um arranged in parallel and 
connected together at one end (Fig. Ufa)). The two reser- 
voirs are separated by a smaller channel of width 10 fxm 
and length 20 /um. Electrodes were fabricated beneath 
the reservoir microchannels and are denoted as the left 
and right reservoir electrodes respectively. A split-gate 
electrode was fabricated at the base of the small central 
channel (Fig. EJb)). The geometry of the central channel 
is shown in Fig. [He); the split gate was 1 /im long and 
separated by a gap of 2.8 /um. The 400 nm gap between 
left and right reservoir electrodes was placed 10 /um to 
the right of the split gate in order to avoid the distortion 
of the potential profile of the constriction region. A 1.5 
/im thick layer of hard baked photoresist, which defined 
the microchannel depth d, separated the lower electrodes 
from the guard electrode. All metal layers were made of 
gold (65 nm thick) on top of a thin (15 nm) titanium 
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FIG. 3. (Color online) Finite element modeling results, (a) 
Contour plot of the calculated electrostatic potential in the 



V st = 0V and V r = +1.0 V. The 



central channel for V g , 
darker areas correspond to regions of lower potential energy 
for electrons, (b) Electrostatic potential in the x direction 
along the channel (y = (jm) and the y direction across the 
channel (x = ^m), for V gu = V, V r = 0.3 V and Vgt =+1.3 
V (red), +0.3 V (blue), -0.7 V (green), -2.3 V (dark yellow). 



layer which was deposited in order to promote adhesion. 

The sample was placed in an experimental cell ap- 
proximately 0.5 mm above the bulk surface of super- 
fluid 4 He at 1.25 K. The potentials V gu , V r and V gt 
were applied to the guard, reservoir and split-gate elec- 



+1.0 V, the surface of the helium was 



trodes, respectively. Under the bias conditions V g ,, 

V, V r = V gt 

charged by thermionic emission from a small tungsten 
filament placed a few mm above the sample. A small AC 
voltage Vi n of frequency 200 kHz was superimposed on 
the right reservoir electrode in order to drive electrons 
between the two reservoirs, through the central channel. 
The current / and conductance of the electron system 
G were determined by making a phase-sensitive mea- 
surement of the voltage capacitively induced on the left 
reservoir electrode, in reference to the standard lumped- 
circuit model^i. 

To aid the understanding of the electrostatic potential 
profile of the device, a finite element model of the central 
channel was developed^. The results of the modeling 
are shown in Fig. |3][a) for V gu — V gt = V, V r = +0.3 
V. When V g t is more negative than V r , a saddle-point 
potential is formed on the helium surface at the point 
(x = 0, y = 0), with a maximum in potential for electrons 
in the x direction along the channel and a minimum in 
the lateral y direction. Fig. [3jb) shows how the poten- 
tial profile develops as V gt is changed. For V gt =-2.3 V 
the potential at the center of the saddle-point is more 
negative than the guard electrode potential. We assume 
that the electrostatic potential energy of the electron sys- 
tem may not exceed V gu , as this would lead to electrons 
escaping from the reservoirs, onto the thin helium film 
covering the guard electrode, through which they should 
rapidly drain away^. Therefore, when the potential at 
the center of the saddle-point is more negative than the 
guard electrode potential, the resulting potential barrier 
between the two reservoirs should block electron trans- 



port through the constriction. 



III. RESULTS AND DISCUSSION 

Before discussing experimental results, we first develop 
an electrostatic model of the device, following the ap- 
proach presented previously^. Note that in this model, a 
more positive voltage corresponds to a lower potential en- 
ergy for an electron. We estimate that the change in the 
depth of the helium at the center of the 20 /mi-wide reser- 
voir microchannels should be less than 0.1 /j,m. Therefore 
we do not take the curvature of the helium surface into 
account here. We begin by assuming that the electron 
system may be considered as a charge continuum. The 
electrostatic potential energy of the electron system in 
the reservoirs, V e , depends on the reservoir electrode volt- 
age V r and the electron density. Therefore we may write 
V e = —en s d/eeo + V r . The results of the finite element 
modeling analysis (Fig. [3]) indicate that a saddle-point 
potential is formed at the constriction between the two 
reservoirs when V g t is more negative than V r . At ap- 
propriate values of V g t a potential barrier between the 
two reservoirs may be formed. The saddle-point po- 
tential may be written as V(x,y) = Vj, + ^ax 2 — \by 2 
where V> is the potential at the center of the saddle- 
point and a and b are constants. We consider that for 
— eV e > —eV> electrons may pass over the potential bar- 
rier, allowing transport through the constriction, whereas 
for — eV e < —eVb electron transport is blocked. The con- 
dition —eV e — —eVb therefore defines the threshold of 
current flow through the constriction. This condition is 
depicted schematically in Fig. EJa). 

The influence of the reservoir, split-gate and guard 
electrodes on the potential barrier may be estimated by 
considering the region of space at the center of the con- 
striction to have some capacitance to the three electrodes, 
C r , C g t and C gu , as well as a stray capacitance to the sur- 
roundings C s . The total capacitance of the region is then 
described by Cs = C r + C gt + C gu + C s and we can define 
constants to represent the relative strength of coupling 
from each electrode to the barrier region as a — C r /Cs, 
/3 = C gt /C-z, 7 = C gu /CYi and a = C s /C^. Experimen- 
tally it was found that, on increasing V r , V g t and V gu by 
+100 mV simultaneously the threshold of current flow 
also increased by exactly +100 mV (data not shown), in- 
dicating that C s was in fact negligibly small (a — 0) and 
that a+/?+7 = 1. The potential at the center of the con- 
striction may then be written as V> = caV r + (3V g t + "fV gu . 

We now consider the case in which the reservoir volt- 
age is set progressively more negative with a fixed bias 
applied to the split gate and guard electrode (here we as- 
sume Vi n to be small). This is depicted schematically in 
Fig. |UJa-c); from top to bottom V r goes from positive to 
negative bias. In Fig. |4ja) we see that, at a certain posi- 
tive reservoir bias, the height of the potential barrier just 
meets the potential energy of the electron system and the 
current is 'pinched-off '. Making V r more negative (Fig. 
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FIG. 4. (Color online) Two electrostatic models of the system, 
(a) Schematic diagram of the electrostatic potential across 
the device through the central channel, V{x, 0), for the case 
where V e = 14, assuming that the electron system forms a 
charge continuum, (b) As V r becomes more negative, V e be- 
comes negative more quickly than Vb and V g t must be set more 
negative to 'pinch-off' the current, (c) Eventually V e = V gu 
and when V r becomes more negative electrons are lost to the 
guard electrode. Because V e now remains constant, V g t must 
be set more positive to maintain the threshold condition, (d) 
Schematic representation of the constriction under a granu- 
lar charge model in which the constriction is modeled as a 
small square of the helium surface of area A — w\. The num- 
ber of electrons in the constriction N c — n s wl. Under this 
model, the current passing through the constriction should be 
suppressed when N c < 1. 



0|b)) causes both V e and Vb to become more negative. 
As V e remains below the guard potential electrons do 
not escape from the channels onto the thin helium film 
above the guard; there is no change in the electron den- 
sity and V e simply changes by the change in V r . However, 
the change in the barrier height is smaller due its addi- 
tional dependence on the split-gate and guard electrodes 
(a < 1). Therefore a negative change in V r causes Vb 
to become more positive relative to V e and current flows 
across the barrier. By setting V r more negative still (Fig. 
|4{c)) electrons are eventually lost to the guard leading to 
a reduction in n s until V e is equal to the guard potential. 
Now as V r becomes more negative, V& becomes more neg- 
ative with respect to V e and the current may eventually 
be suppressed once more. 

We denote the value of V g t for which — eV e — —eVb as 
V gt . Following the model described above, expressions 



may now be derived to describe the dependence of V gt 
on V r , for constant V gu . From the threshold condition 
— eV e = —eVb, for the case where — eV e < ~eV qu we have 



ee 



+ V r = aV r + /3V g ° t + jV gu , 



and for the case where — eV e 



V a , 



aV r + (3V Q ° t + 7 V 6 



Rearranging (1) and (2) gives respectively, 
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Because the separation between the split-gate elec- 
trodes (2.8 /xm) is comparable to the inter-electron spac- 
ing (~ 0.3 /xm for n s = lx 10 9 cm -2 ), we expect only a 
small number of electrons to be in the constriction region 
when the system is close to the threshold condition. To 
estimate this number, we model the center of the con- 



striction as a small square of area A 



where 



is the effective width of the constriction. As the FEM 
calculation indicates that the constant a is small com- 
pared to b, we will assume that the potential in the x 
direction is flat, whilst the electrons are confined in the 
y direction by a parabolic potential V(y) = Vb — ^by 2 , as 
pictured schematically in Fig. HJd). The maximum lat- 
eral displacement for electrons in the constriction, y maX i 
satisfies the expression V(y max ) = Vb - \by 2 max = V e . 
Substituting y max = w c /2 gives 



-(Vb ~ V e ) 



(5) 



Note that for the threshold condition — eV e — — eVb, w c = 
0. Assuming a parallel-plate capacitor approximation, 
the charge density in the area A may be written as 



ed 



(V(y)-V e ) , 



and the total number of electrons in A is therefore 



££o 
ed 



(V(y) - V e ) dxdy 



(6) 



(7) 



After integration over the limits as given by Eq. ([5]), we 
obtain the result 



(8) 



Although for the electron liquid electrons are not local- 
ized as in the Wigner crystal, we assume that electrons 
are distributed evenly over A, effectively forming a series 
of rows across the constriction. The number of electrons 
lying in the y direction across the constriction may then 
be estimated as N y = y/N c . For constant V r and V gu , we 
also have the relation Vb — V e = P(V g t — V gt ). This gives 
an expression relating N y to V g t as 



Ny = 



2 8ee 

3 bed 
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V g ° t ) 



(9) 



The number of electrons across the constriction should 
therefore increase linearly with increasing V g t above the 
current threshold. 

We suggest that the granularity of charge may cause 
some deviation from the behavior of the system as ex- 
pected under the charge continuum model. From Eq. 
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©, the value of V g t for which one row of electrons may 
be formed in the constriction (N y = 1) is 



v 9 \ 



gt 



3 bed 
2 8e^ 



(10) 



For Vg t < Vgt < Vg t the number of electron rows across 
the constriction is less than 1. Therefore, as no electrons 
are present in the constriction, current flow should be 
suppressed when V gt = V gt rather than at the threshold 
condition assumed under the charge continuum model 
V g t = Vg t . We may now correct Eqs. $5$ and (@| to give 
Vg\ for changing V r at constant V gu . For the case where 
-eV e < 



-eV gu we have 



1 3 b(V r )ed 
P\2 8ee 
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+ 7Vg U 



, (11) 
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and for the case where — eV e 



-eV„, 
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Note that the value of b, which describes the parabolic 
lateral confinement at the constriction, is dependent on 
the potential applied to the reservoir electrodes (assum- 
ing constant V gu ). The FEM calculation of the potential 
profile shows this dependency to be linear. The correc- 
tive terms therefore introduce a non-linear dependence of 
V g \ on V r . However, in the following discussion, we will 
assume that, for constant V r , small changes in V g t cause 
a negligibly small change in b. 

We have performed measurements in order to deter- 
mine whether the charge continuum model or the granu- 
lar charge model best describes the dynamics of electrons 
in the device. We denote the experimentally determined 
value of the split-gate voltage for which the current is 
suppressed as Vffi ■ Figure 02a) shows the dependence 
of the current / flowing through the central channel as 
Vgt is varied for different values of the driving voltage 
Vin- The current decreases as V gt is swept negative and 
finally is completely suppressed at the threshold voltage 
V'/ 1 , which for all cases is more negative than the reser- 



-1 V. The linear dependence of V*/ 1 



"gt 

voir voltage V r — -,- ± v . j. mnoi uc F wiu^ ui v gt 
on Vi n is shown in Fig. [51(b). The threshold for current 
flow occurs at more negative split-gate bias as the driving 
voltage is increased. 

Measurements of SSE current flow over a potential bar- 
rier formed by a split-gate electrode in a similar device 
have already been reported^. There, for values of Vi n 
up to 100 mVpp, it was shown that current flow across 
a potential barrier could be induced by increasing Vt n 
until the potential energy of the electron system could 
overcome the barrier, allowing electrons to be transferred 
between the reservoirs. Here we observe the same behav- 
ior. The dependence of Vffi on Vi n is discussed in more 
detail later in this section. 



FIG. 5. (Color online) (a) Peak current I measured as a 
function of V g t for different driving voltage amplitude Vm- 
(b) The current threshold Vgt against Vin- (c) Conductance 
of the electron system G as a function of V g t for different 
Vin- Above 4 mVpp each data set is shifted vertically for 
clarity, (d) Constriction conductance G c for V% n = 4 mV pp . 
G c increases in a series of smoothed steps. The grey dotted 
line is a guide to the eye indicating values of V g t at which, 
under the granular charge model, the number of electrons 
across the constriction increases by 1. The threshold V g f lies 
between N y = and N y — 1. The solid grey line is a guide to 

3 

the eye of functional form G c ~ V g \ starting at V g t = 0.25 V. 



In Fig. 03c) we show the corresponding conductance G 
for each current measurement. Above V g t ', G increases 
sharply for all values of Vi n . Considering the saddle- 
point potential shown in the modeling results in Fig. [31 
both the depth of the potential at the center of the con- 
striction (which is proportional to the electron density 
in the constriction under the charge continuum model), 
and the effective width w c of the constriction, should in- 
crease as V g t is swept positive. Both effects should lead 
to an increase of conductance at the constriction. For 
high values of V g t , where the split-gate electrode is more 
positive than the reservoir electrode and no potential bar- 
rier is expected to exist, G rises further, presumably due 
to the continuing increase of the electron density in the 
constriction region. 

For still higher values of V g t, the conductance even- 
tually saturates. For V m = 4 mV pp , G = 2.30 Mft -1 
at Vgt = 3.0 V. As the total conductance of the elec- 
tron system no longer depends on the split-gate voltage 
in this region, we assume that the resistance of the con- 
striction region is small compared to that of the electron 
system in the reservoirs. Therefore, the saturated con- 
ductance can be attributed to that of the reservoir re- 
gion, and by subtracting the value of the resistance at 
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this split-gate voltage, R = 1/G = 0.434 Mil, from the 
values of the resistance over the entire split-gate sweep, 
an approximation of the constriction resistance R c can be 
calculated. The corresponding conductance of the con- 
striction, G c = 1/ R c , is shown in Fig. [5jd), for V in = 4 
mVpp. We see that above the conductance threshold, 
G c increases in a series of steps. The steps are not a 
series of sharp rises and flat plateaus; rather, they ap- 
pear smoothed-out. We add to the plot a guide to the 
eye (grey dotted line) indicating the manner in which a 
sharp step pattern fits the data. We align each sharp 
step with the maximum in the gradient at each increase 
in G c . The spacing between the steps is AV gt = 250 mV. 

We suggest that the increasing number of electron rows 
across the constriction could lead to steplike increases 
in the constriction conductance. Such behavior is ob- 
served in other classical many-body systems with long- 
range interactions, such as pedestrians moving through 
bottlenecks^.. Here, such an effect is essentially the re- 
sult of Coulomb blockade at a single constriction; for 
the case in which N y = 1 and one electron row occupies 
the constriction, Coulomb repulsion prevents other elec- 
trons from passing through and electrons may only pass 
through in a single row. As the constriction is opened 
to the point where N y — 2, the additional conduction 
channel should cause an increase in the constriction con- 
ductance as electrons may now pass freely through the 
constriction, side-by-side. The steplike increase in con- 
ductance resembles the quantized conductance steps ob- 
served in quantum point contact devices^. From the 
FEM calculation results shown in Fig. [3Jb), we esti- 
mate that the spacing of the energy subbands for lateral 
motion at the constriction to be ~ 0.1 meV, which, as 
shown below, is much smaller than the change in Vb as- 
sociated with each conductance step, and is also smaller 
than Vi n . We therefore conclude that the origin of the 
steps observed here, due to the Coulomb interaction be- 
tween electrons at the constriction, is quite different to 
the case of the quantum point contact, of which our de- 
vice acts as a classical analogue. 

The smoothing of the steps in G c indicates that nei- 
ther the charge continuum model, nor the granular charge 
model, describes the system precisely. We assume that 
the sharp steps fitted to the data in Fig. EJd) correspond 
to the values of V g t where, under the granular charge 
model, the number of electrons able to pass simultane- 
ously through the constriction increases by 1 . By extrap- 
olation of these fitted steps, we can obtain an estimate 
of Vg t = 0.25 V. As described above, the constriction 
conductance should increase with increasing w c and as 
\b becomes more positive. The FEM calculation shows 
that w c ~ \jV g t whereas Vb ~ V gt - We therefore naively 
expect the conductance of the constriction to vary as 

3 

G c ~ V„t , under the charge continuum model. In Fig. 

EKd) we plot the function G c = 1.04 x 10 7 (V gt - V%)* 
(solid grey line). We see that this function describes 
G c reasonably well, other than at points along the curve 
where the conductance is suppressed, which we attribute 



to deviations from the continuum model due to the gran- 
ular nature of charge. The agreement becomes closer for 
higher values of Vgt as the step features are lost, pre- 
sumably as the number of electron rows increases and 
the electron system at the constriction better approxi- 
mates a 2D charge continuum. We therefore conclude 
that whilst the granularity of charge causes observable 
deviations from the charge continuum model, and causes 
the current to be suppressed above the expected thresh- 
old Vg t , the simplified approach depicted in Fig. 2(d) 
does not describe the system accurately. Correspond- 
ingly, the experimentally observed threshold of current 
flow, V*t — 0.456 V, lies between the values which we 
have estimated correspond to V gt and V gt . 

We suggest that our experimental observation of the 
step-like increase in the conductance of the constriction 
may be verified by molecular dynamics simulations of 
classical charge systems at potential bottlenecks^. In- 
deed, recent simulations of a system very similar to ours 
reproduce the smooth steplike increase in G c , and show 
that each step is related to an increase in the number 
of electron rows across the constriction^!. These sim- 
ulations also indicate that temporal fluctuations in the 
potential of electrons at the constriction, which are of 
thermal origin but essentially due to electron-electron in- 
teractions, and are not considered in our mean-field ap- 
proach, cause the smoothing of the conductance steps. 
These fluctuations also cause V g t to lie at vales of V gt 
more negative than V gt as, even for the case N y < 1, there 
still exists some probability that the potential of electrons 
may be raised in order to overcome the barrier at the con- 
striction, thus allowing transport. For our experiments, 
the driving voltage V% n causes an additional modulation 
of the electron density at the constriction over each AC 
cycle. In Fig. [5^c) the smoothing of the step-like fea- 
tures increases with increasing V% n , as discussed in more 
detail later in this section. However, we note here that 
extrapolation of the data shown in Fig. [5jb) to the limit 
Vi n — mVpp yields a value of V g t ~ 0.47 V, which is a 
small change from the value recorded for Vi n = 4 mV pp 
and is still more negative than the value of V g \ = 0.5 
V. Further comparison of our experimental results with 
numerical simulations is required to fully understand the 
dynamics of the electron system close to the conductance 
threshold. 

We consider that the modulation of the potential of 
the electron system due to the applied driving voltage 
may cause additional non-linear behavior to appear in the 
electron transport properties close to the current thresh- 
old. To investigate such effects, the magnitude of the 
second harmonic component of the AC voltage induced 
on the left reservoir electrode, i?2/, was measured. In the 
case that the SSE AC current flow is perfectly sinusoidal 
in response to the driving voltage, the second harmonic 
component is zero. However, if the conductance should 
vary over each AC cycle, the current signal should be- 
come distorted causing higher frequency components to 
appear. In Fig. UJa) we see that i?2/ indeed rises as V g t 



FIG. 6. (Color online) (a) Magnitude of the second harmonic 
component of the voltage signal, f?2/ , as a function of V g t for 
different driving voltages, (b) Processes giving rise to distor- 
tion in the AC current through a saddle-point potential. As 
the potential across the constriction V y is parabolic, the ef- 
fective width of the constriction w c oscillates over each AC 
cycle of Vin, causing the conductance of the constriction re- 
gion to vary over each cycle, (c) Close to the threshold of 
current flow, electrons may only pass intermittently across 
the potential barrier formed at the constriction. 



is swept negative and reaches a maximum close to V g f. 
The magnitude of i?2/ increases with the driving voltage 
V m . 

In Fig. |6jb) and (c) we depict schematically two pro- 
cesses which may give rise to the distortion of the electron 
current through a saddle-point potential as the electron 
energy V e is modulated. In the first case (Fig. [^b)), 
when V e is close to the bottom of the parabolic poten- 
tial V y , the effective width of the conductive channel w c 
changes over the AC cycle which should lead to a time 
dependence of the conductance G. In addition, the effec- 
tive depth of the potential for electrons at the constric- 
tion also varies, which should lead to a variation of the 
electron density, and so conductance, over each cycle. In 
the second case (Fig. EJc)), as the potential maximum 
in V x is raised and becomes higher than V e — Vi n /2, elec- 
trons are expected to flow only intermittently across the 
barrier. 

In both of these cases, the degree of distortion in the 
current should reach a maximum close to the threshold of 
current flow. We therefore find the observed increase in 
i?2/ close to the current threshold to be consistent with 
transport through a saddle-point potential. We also note 
that i?2/ appears to rise in a series of weak peaks or steps 
which appear to mirror the step-like decrease in conduc- 
tance observed in Fig. G3Jc). We again find this to be con- 
sistent with our model. For values of V g t where dI/dV g t 
is large, modulating V e gives rise to a large change in 
conductance over each AC cycle, and so the value of 
should be large. Each step in conductance should there- 
fore be accompanied by a peak in i?2/- Indeed, the step- 
like feature in the conductance at V gt = 0.75 V (Fig. 
[He)) is accompanied by a weak second peak in i?2/ at 
the same split-gate voltage. However, the non-linear re- 
sponse of the SSE system to a driving field has been pre- 
dicted by Saitob~, and observed experimentally 39 . Wc 
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FIG. 7. (Color online) Measured values of for different 
values of Vr, which was varied from positive to negative val- 
ues, and Vi„ — 8 mV pp . The measurement was performed 
five times, as indicated by the different colors. The helium 
surface was charged at the start of each measurement. The 
dotted and dashed lines correspond to Eqs. ([3jl and (j4| re- 
spectively. As explained in the text, for the data sets 1-5, 
n s = 0.70, 0.50, 0.99, 1.28, 1.48 x 10 9 cm" 2 respectively. 

therefore note that the distortion in the current signal 
may also be due in some part to the intrinsic non-linear 
transport properties of the SSE system rather than solely 
due to the geometry or potential profile of the device. 

The dependence of V g ^ on V T was investigated by mea- 
suring the current threshold for decreasing values of V r , 
from 1.5 V to V in 50 mV steps, for V in = 8 mV pp . The 
results of 5 such measurements, taken on different days, 
are shown in Fig. [71 In all 5 cases the split-gate voltage 
required to suppress the current initially becomes more 
negative as the reservoir potential is made more nega- 
tive. Then, at a certain reservoir potential, the trend is 
reversed; the threshold moves to more positive values as 
the reservoir potential becomes more negative. In both 
cases the relationship between V^* t and V r is approxi- 
mately linear. These results are found to be in agreement 
with the behavior predicted by our electrostatic analysis. 
Starting at highly positive reservoir electrode bias, V g t 
must initially be made more negative to suppress the 
current flow each time V r is set more negative, to com- 
pensate for the reduction in the barrier height relative to 
V e - Then, as electrons are lost to the guard, n s decreases 
and V e remains constant, the split-gate voltage must be 
made more positive to allow current flow with each step 
in V r , due to the relative increase of the barrier height. 
The intersection between the two linear regions in the 
data marks the point at which electrons begin to escape 
to the guard. This point is different for each data set in- 
dicating that the initial surface density in each case was 
different. The intersection was not observed to occur for 
values of V r greater than +1.0 V, the voltage at which 
the electron reservoirs were charged. 

By making linear fits to the data presented in Fig. 
[7j we may derive values for a, ft and 7 based on the 
charge continuum model, using Eqs. ((3]) and Q, and 
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the relation a + /3 + 7 = 1. The values of the coupling 
constants were determined for each of the data sets. The 
similarity in the gradient of each data set indicates the 
stability of the coupling constants, the average values of 
which are shown in Tabic [I] 

We may also estimate a, (3 and 7 taking into ac- 
count the fact that the experimentally determined cur- 
rent threshold does not correspond exactly to V gt . The 
results in Fig. [5] show that V g t appears to lie closer to 
Vg t , the threshold under the granular charge model, as 
given by Eqs. (fTTj) and (fl"2"|) . Whilst the constant b, 
which determines the deviation of V gt from V g \, may be 
estimated from the FEM modeling of the device, here we 
adopt a more straight-forward approach to take into ac- 
count the deviation of Vffi from V gt . From Fig. [5] for 
V r = 1.0 V, Vgf - V° t = 0.16 V. As discussed in more 
detail below, the reservoir electrode potential is equal to 
the guard electrode potential for V gu 0.5 V, due to an 
offset in the potential of the guard electrode. Because 
no parabolic confinement exists when V r — V gu , for this 
condition we expect 6 = and therefore Vgt — V gt — 
V. As the correction to V gt should display a square-root 
dependence on V r , the value of V gt can be estimated for 
each value of V r using the expression 

V g ° t - V g f - 0.16^-^ . (13) 

The average values for a, /3 and 7 calculated using the 
values for V gt given by Eq. (fT3"|) are also listed in Table 
III 

The coupling constant j3 may also be estimated from 
dependence of Vgt on V in shown in Fig. [SJb). We as- 
sume that, on increasing Vi n , the increase in the bar- 
rier height required to suppress the current is equal to 
the corresponding increase in the maximum energy of 
the electron system in the right reservoir, Vj n /2. Be- 
cause, for fixed bias on the reservoir and guard electrodes, 
AVfc = flAVgt we may derive /3 from the gradient in 
Fig. Mb) as dVgt/dVin — —1/2(3. Similar measurements 
were made by sweeping the reservoir and guard potentials 
(data not shown) in order to determine the dependence 
of the thresholds V r th and V g f on Vin ■ Under the charge 
continuum model, the constants a and 7 may be given 
by dV r th /dV m = 1/2(1 - a) and dV g f/dV m = -l/2 7 . 
The results of these measurements (denoted dV /dVi n ) 
are shown in Table U as well as the corresponding values 
calculated from the finite element model. 

From the results shown in Table [H we see that, com- 
pared to the values calculated using Eqs. ©, (0]), the 
change in a, /3 and 7 when including the correction given 
by Eq. (fT3"|) is small. These values are also close to those 
given by measuring the dependence of V th on Vi n . All the 
experimentally determined values are in relatively good 
agreement with the FEM calculation. For all four results, 
a is largest and so the reservoir electrode dominates in 
determining the height of the potential barrier. This con- 
firms that in our experiment, close to the conductance 



threshold, the electrons are indeed passing through the 
region between the split-gate electrodes and above the 
reservoir electrode. 

The determination of the constant /? allows the ex- 
pected change in V g t to add one electron row across the 
constriction to be calculated using Eq. (ITU1) . We use the 
FEM model to obtain an estimate of b — 2.5 x 10 11 Vm~ 2 
for V r = 1.0 V and ^=0.62 V. Using /3 = 0.16 we find 
that Vg\ — Vgt = 224 mV. This value is in good agree- 
ment with the experimentally observed step separation, 
AVgt = 250 mV for V r = 1.0 V (Fig. EJd)), indicat- 
ing that the step-like increases in G c are indeed due to 
increases in the number of electrons across the constric- 
tion. The change in split-gate voltage AV gt is related to 
the corresponding change in V\, by the constant /3. For 
f3 = 0.16, AT4 « 40 mV, which is comparable to V in . We 
therefore conclude that the increase of smoothing of the 
steplike features with increasing Vi n may be due to the 
modulation of the electron density at the constriction, 
which causes transport features arising from the discreet 
number of electrons across the constriction to be lost. 

The data shown in Fig. [7J was taken with a voltage 
of V applied to the guard electrode. From Eq. @ it 
is then to be expected that V g f = V for V r = V. 
However, the intercept on the V g f axis for each data set 
is typically +3.6 V indicating a true value of the guard 
potential of V gu ~ +0.62 V. This apparent offset is con- 
firmed experimentally for each of the data sets; below 
a value of V r s» 0.5 V the electron signal was lost com- 
pletely indicating that at this point V r ps V gu . The offset 
is seen to vary over a small range; for data sets with a 
higher intercept on the V g f axis the electrons were lost 
at more positive values of V r . Despite thermally cycling 
the device to room temperature and making efforts to 
ensure correct grounding of the experimental wiring the 
offset remained. However, the stability of the offset over 
many hours or days allowed consistent measurements to 
be made. Voltage offsets in mesoscopic devices may be 
caused by contact potential differences, thermoelectric ef- 
fects or surface charging effects^,. The cause of the offset 
observed in this experiment is not yet clear. Taking the 
offset into account, the estimated initial surface electron 
densities for the data sets 1-5 shown in Fig. [7J are 0.70, 
0.50, 0.99, 1.28, 1.48 x 10 9 cm" 2 respectively. 

The melting temperature of the 2D Wigner solid de- 
pends on the electron density as T m = 0.225 x lQ & ^/nJ&. 
For the case of a potential offset on the guard potential 
of +0.62 V, the saturated electron density achieved when 
charging the reservoirs with V r = +1 V is n s — 1.48 x 10 9 
cm~ 2 which gives a value of T m — 0.866 K. In Fig. [5] we 
show the temperature dependence of the resistance R of 
the electron system for V r — V g t — +1.0 V. At 1.23 
K the current threshold was measured (inset). In refer- 
ence to Fig. [7J the highly negative value of the threshold 
Vgl L — —0.29 V indicates that the density was close to 
saturation. The increase in R below 1 K is attributed to 
the formation of the 2D Wigner solid. The localization 
of electrons in the Wigner lattice leads to the formation 
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FIG. 8. (Color online) Resistance against temperature for the 
case in which the constriction is widely opened (V g t = +1 V) 
for Vin = 5 mVpp . Inset: / as a function of V g t for the same 
electron density, measured at 1.23 K. 

of a small depression, or dimple, in the helium surface 
beneath each electron, which increases the electron effec- 
tive mass, and so resistivity of the system^. Recent ex- 
periments have investigated the transport of the Wigner 
solid in a microchannel geometry where the decoupling 
of the electron lattice from the dimple lattice at high 
driving fields leads to a highly non-linear response^. In 
our measurement, a sharp increase in the resistance is 
observed at T s» 0.875 K which may correspond to the 
melting temperature of the electron system, in agreement 
with the predicted value of T m for the case in which the 
guard potential is offset. 



constriction formed by a split-gate electrode. The elec- 
tron current may be suppressed by sweeping the voltage 
of the split gate negative. The threshold for current flow 
was dependent on the DC voltages of all the device elec- 
trodes and the AC driving voltage applied to the electron 
system, as well as the electron density. Step-like increases 
in the conductance of the electron system as the split- 
gate voltage was swept positive were found to be due to 
increases in the number of electrons able to pass simul- 
taneously through the constriction. The device therefore 
acts as a classical analogue of the quantum point contact. 
Our results are in good agreement with a simple model of 
the device developed with the aid of finite element anal- 
ysis software in which a saddle-point potential profile is 
formed at the constriction. Comparison with this model 
reveals that a potential offset on the guard electrode of 
the device plays a crucial role in determining the poten- 
tial profile of the sample. Such detailed characterization 
of microfabricated samples for electrons on the surface of 
liquid helium is an important step towards the realiza- 
tion of more advanced mesoscopic devices such as single 
electron devices and quasi-one dimensional wires. 
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TABLE I. Coupling constants a, ft and 7 as calculated by the finite element analysis and the values measured experimentally 
from the results shown in Fig. [5] and Fig. [7] 



Coupling constant 


Calculated (FEM) 


Measured (Eqs. ©, ©) 


Measured (Eq. ([T3|)) 


Measured (dV th /dV in ) 


a 


0.75 


0.77 


0.79 


0.83 


P 


0.10 


0.16 


0.15 


0.10 


7 


0.15 


0.07 


0.06 


0.07 



